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We investigate the temporal Wilson loop using the Hamiltonian approach to Yang-Mills theory. 
In simple cases such as the Abelian theory or the non-Abelian theory in (1+1) dimensions, the 
known results can be derived using unitary transformations to take care of time evolution. Al- 
ternatively, the exact solution can also be found in Coulomb gauge using the exact ground state 
wave functional which is known explicitly in these simple cases. The Coulomb gauge technique 
can also be applied to the more reaUstic case of Yang -Mills theory in (3+1) dimensions, where 
one has to rely on the approximate vacuum wave functional obtained e.g. in recent variational 
approaches. We use this formulation to compute the temporal Wilson loop in (3+1) dimensional 
Yang-Mills theory, and find that the Wilson and Coulomb string tension agree within this approx- 
imation scheme. Possible improvements of these findings are briefly discussed. 
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1. Introduction 

Recent efforts to study Yang-Mills (YM) theory are largely based on Coulomb gauge, be- 
cause the Gribov-Zwanziger picture of confinement [|l|, ^ becomes particularly transparent in this 
framework: the only constraint on the physical states in a Hamiltonian approach based on Coulomb 
gauge is Gaufi' law, which in turn can be resolved exactly so that any normalizable wave functional 
can be used as a physical state. 

Such a fonnulation naturally lends itself to variational methods, which are usually based on 
(modified) Gaussian Ansdtze for the vacuum wave functional, see e.g. ^ and references therein. 
The solutions agree with the Gribov-Zwanziger scenario and, when combined with the horizon 
condition, lead to a coherent picture of strongly coupled YM theory, which is mainly dominated by 
a diverging gluon energy in the infrared, a lineai^ly rising Coulomb potential between static quarks, 
and a strongly enhanced Faddeev-Popov ghost propagator in the infrared. The latter condition 
implies that the YM vacuum behaves as a perfect colour dia-electric medium [^], i.e. a dual su- 
perconductor. All these findings agree qualitatively with the results from recent lattice simulations 
[^], although there are still some discrepancies in the quantitative details, in particular for the ghost 
sector [^. 

However, this simple picture is subject to the caveat that the Coulomb potential is only an 
upper bound for the tme potential between heavy quarks [^. For a complete description, one 
has to address the physical potential from large Wilson loops directly. In the continuum, this is 
complicated by path ordering and large self-energy contributions which obscure the extraction of 
the true potential from large Wilson loops [^]. 

In this talk, I report on a recent approach [ p^ to compute the Wilson potential in the Hamil- 
tonian approach in Coulomb gauge. A careful analysis of GauB' law and the physical degrees 
of freedom allows for an elegant way to isolate and remove the self-energy contributions to the 
Wilson loop, while unitary transformations are applied to take care of time evolution. As specific 
examples, we study quantum electrodynamics (or rather Maxwell's theory) in D = (3 + 1) space- 
time dimensions, and YM theory in both D = (1 + 1) and D = (3 + 1). In the latter case, only an 
approximate solution can be found and I briefly discuss possible ways to improve our findings. 

2. The Wilson loop in the Hamiltonian approach 

We study SU {N) YM theory in euclidean space using the Hamiltonian approach in Weyl 
gauge, Aq = 0. The physical energy states of this formulation are stationary wave functionals 
*F„[A(x)] • e^'^"' which obey the Schrodinger equation based on the Hamiltonian 

H = ^ Jdxn^{x) + ^ JdxB^. (2.1) 

Here, A"(x) is the gauge connection, n"(x) = —id/5A"{x) its conjugate momentum, g the (bare) 
coupling constant and B" = V x A" — ^/"^'^A* x A"^ the non-Abelian magnetic field. 

To proceed, we must fix the residual time-independent gauge symmetry left after imposing 
Weyl gauge, for which we choose the Coulomb condition, V • A = 0. This adds gauge-fixing terms 
to the Hamiltonian (which we discuss below) and imposes Gaufi' law as a constraint on the physical 



2 



Temporal Wilson loop in the Hamiltonian approach in Coulomb gauge 



Markus Quandt 



states: In the sector with n static chai^ges located at positions xi . . . ,x„, we have n additional colour 
indices on the wave functionals and GauB' law reads (with generators in the right representation) 

(z) ,..,■„ (A ; XI , . . . , x„) = - I £ I 7^-^,, 5 (z - X,) n 5,„,.,„ I (A ; xi ,.. . x„) . (2.2) 

U=0{k,} m^e ) 

The complete physical Hilbert space thus decomposes into orthogonal charge sectors characterised 
by the data from the external charges. 

Next we consider a rectangular temporal Wilson loop of extension {T x R). Since Aq = 0, we 
have W = tr{ U [A] {R,0;Ty -U [A] (0, 0; 0) } , where U [A] (y, x; f ) is the parallel transporter x y at 
fixed time t, in the background of the gauge potential A. Using euclidean time evolution and the 
abbreviation U (x) = ?7[A](x,0;0), we find 

(W) = tr (0 1 U(Ry e-^"-'^"'^^ U{R) 1 0) = tr (R | £-^^(^-^0) | _ (2.3) 

The salient point is now that the vacuum is gauge-invariant, while the Wilson state \K) =U (R) |0) 
is in the ^^-sector according to eq. (2^). Inserting a complete set of energy eigenstates from all 
sectors, only the <7^-sector will thus contribute and we have 

(W)=tr£(0|[/HR)|«),^-e-(^-^«)^-,^(«|[/(R)|0)=££'-^(^""-^«)|KR (2.4) 

n n 

where = tr(n^n). In the limit of large Euclidean time extensions, we project out the Wil- 
son potential as the difference between the ground state energies with and without static quarks, 
Ef{R)-Eo = -limr^„ T'^ ln{W) . 

The next step is to remove the Wilson line from the state |R) by means of a unitary transforma- 
tion based on the parallel transporter U (R). Under such a transformation, the conjugate momentum 
(and thus the Hamiltonian eq. (2.1)) acquires an additional term, n"(x) — )• n"(x) = n"(x) + £?(x), 
where 

e"ix)= rjy5(x-y)?7(y)t/r[/(y). (2.5) 

JO 

As a consequence, the Wilson loop is now computable as a zero-charge vacuum expectation value, 

{W)=tr{0\e-^^"-^'>'>\0) (2.6) 

where H differs from eq. (2J_) by the momentum shift ^ 

To illustrate the problems in extracting the true potential, let us consider the contribution to 
the Wilson loop eq. ( ^ ) which is quadratic in the induced electric field £(x), 

2 2 

Vind(/?) = y tr J dx[e%x)f = |-C25f (0) • |R| , (2.7) 

where C2 is the quadratic Casimir operator for the colour group, and the (quadratically) divergent 
prefactor 5j^(0) is due to the infinitely thin Wilson lines. Although Vind is formally confining, it 



This shift complicates the treatment of H considerably, since 
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must be spurious because it exists even for G = U{1). On the other hand, it is also not correct to 
simply drop Vind, because the true potential is hidden underneath the divergence: For G = U (1), 

2 2 

VindiR) = ^tv jdx [£^(x)]2 + |-tr j dx [£||(x)]' , (2.8) 
where the first factor renormalizes the Wilson loop, and the second factor gives the true (Coulomb) 



potential, cf. section 3. 1 



The unitary transformation leading to eq. ( [2.6[ ) only assumes Weyl gauge and thus works with 
or without Coulomb gauge fixing. In cases where the exact gauge-invariant ground state is known, it 
is not necessary to fix the residual gauge freedom since the Wilson loop is gauge invariant. In more 
realistic cases, however, the gauge invariant ground state will not be known and we have to resort 
to approximations based on Coulomb gauge. The required gauge-fixing is virtually impossible to 
effect directly in eq. (^^, except for the Abelian model. A much better strategy is to go back 
to eq. (pj|), perform the standard gauge fixing and resolve GauB' law exactly, which induces in 
H — )• //fix the well-known Coulomb term,^ 

= y / d\x,y) [p"{x)+p%,{x)]/F"\x,y) [p\y)+p%Ay)] , (2-9) 

where ^ is the Faddeev-Popov determinant and F"^ the the so-called Coulomb kernel. The Wilson 
loop now reads 

{W) = tr{R^\e-^^"<'--^'>^R^) (2.10) 

where the new Wilson-like state |Rj^) = U [A±\ (R) |0) is no longer from the qq sector, but rather has 
overlap with the zero-charge vacuum. Any divergences associated with the point-like cross section 
of of the Wilson lines in |Rj^) must therefore be unrelated to external charges: as we will see, 
the extra divergences renormalize the composite operator W — )■ Wr = Z^y(A) • W, but they do not 
affect the physical potential. Whether or not the state |Rj^) in eq. ( 2.10| ) is subsequently removed 



from the Wilson loop by a unitaiy transformation U± (R) as described earlier becomes a matter of 
computational convenience. 

3. Applications 

3.1 Quantum Electrodynamics 

In the Abelian case G = U{1), the exact ground state is known and we can first attempt a 



gauge-invariant treatment based on eq. ( |2.6[ ). The induced electric field 

R 



£(x) = J Jy5(3)(x-y). (3.1) 







is now field-independent, so that the induced potential ( |2.7[ ) could be pulled out of the expectation 
value (W). In view of the discussion above, it is more convenient to first split e = e-^ + and pull 



^As a consequence of the parallel transport, the non- Abelian charges p"(x) = 5(x — R) iT" — 5{x)U (R) iT" U (R) 
are field-dependent. 
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out the longitudinal piece only, (W) = e ^^n^'^) • {0\e ^(^qed £'o)|o) where 

^QED = yd'x [n^(x) + £^(X)] ' + ^ [V X A^(x)] '| . (3.2) 



It is now convenient to reverse the unitary transformation that led to eq. (p77|), but this time with the 
transversal gauge connection A-'- only. As a result, the shift in the electric field is removed from 



the Hamiltonian (|j) and reshuffled into the state \R^) = U[A-^]{R) |0) = U-^{R) |0). Although 



this new state resembles the initial Wilson state, it is gauge-invariant as it depends on A-'- only. 



Inserting a complete set of eigenstates of the standard QED Hamiltonian eq. (3^) (without e^) 
yields 

(W) = e-^^ii • £ 1(0 1 {R)\n) |2 ^-^(^--^o) ^ |(o| [/^ (R) |0) |2 . ^.-^^11 («) , (3.3) 

provided that |R-'-) has non- vanishing overlap with the true vacuum. For Maxwell's theory the 
ground state is Gaussian and the relevant matrix element can be computed exactly. Using a 0(3)- 
invariant UV cutoff A for the momentum integration, we find that formally vanishes in the limit 
A — )■ oo because |R^) has poor overlap with the true ground state [||]. On the other hand, Z^/ is also 
independent of the temporal extension T of the loop and hence cannot contribute to the physical 
potential. 

The conect interpretation is given by the operator product expansion (OPE): Since the loop 
operator W contains products of field operators at arbitrarily close points, we expect short distance 
(A — )• oo) divergences associated with W, on top of the counter terms necessary to render Green 



functions finite. In fact, it has been known for a long time [11, 12] that one overall multiplication 
W ^ Wr = ■ W is sufficient to render (Wr) = e^^^^ finite. From eq. (3.1 1, the true physical 
potential becomes 

VW=.,to4/</.M(x)F = j^-^. (3.4) 

This is just the usual Coulomb potential including the self-energy of the charges. 

The entire derivation could also be repeated in a Coulomb gauge fixed formulation. If we start 



from eq. ( |2.6| ), the gauge fixing gives H — ^ //fix = ^qed + ^| > ^1- ( |3.2| ), and the computation 



becomes identical to the gauge-invariant treatment above. Alternatively, we could also start from 



eq. ( [2.3| ) and resolve GauB' law for the state |R) explicitly, which induces the Coulomb term H — 
^fix = ^qed[A_l] +Hc in the Hamiltonian. This time, the Coulomb potential comes directly from 
He, while the first term in //fix only renormalizes the Wilson loop. 

3.2 Yang-Mills Theory in D = 1 + 1 

This model is (almost) topological and requires a non-contractible spactime manifold to give 
nontrivial results. Since time must be unrestricted in the Hamiltonian formalism, we compactify 
space to an interval [0,L] with periodic boundary conditions and take spacetime to be the cylinder 



M = R X S . For this model, the induced potential eq. (2.7) gives the correct string tension 



^ SU{2) 3 2 

c^i+i = 2"<^2 = g§ ■ (3.5) 
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In view of the above discussion, it is, however, unclear why this should be the correct answer, in 
particular- since Vind = ai+i /? is not periodic or invariant under /? — )• L — /? as expected from the 
compactification. 

To address this issue, we go back to eq. ( |2.3[ ) and fix the minimal Coulomb gauge d\A" = 
and Aj = diag. Taking the colour group G = SU (2) for simplicity, we are thus left with only one 
scalar degree of freedom 



^ = -AiLe[0,n]. 



(3.6) 



The restriction to the compact interval (the fundamental modular region in this case) eliminates 



all residual gauge symmetries. In order to fix eq. (p.3|), we have to determine the gauge fixed 
Hamiltonian in the ^^-sector to which the Wilson state \R) belongs. This has two pieces 

„2 



Hfix = Hf^^ + Hc = //fix + ■ 



(3.7) 



where the non-Abelian charge interaction comes from the resolution of GauB' law. The zero-charge 



Hamiltonian H9 and its spectrum have been studied thoroughly in ref. [13 1. As for the Coulomb 



term, we have to take into account that the external charges p" are field-dependent due to the 
parallel transport. Expanding all colour vectors in a a polar basis (cq = ^3, e± = (ei it 62) / \/2), we 
find 



He 



8^ ^ 

- y ■ 



2fa(0;T?) 



-2ic (R/L)i? 



F^(/?;t?)+cc]} 



(3.8) 



where T(j are the polar Pauli matrices. The Coulomb kernel F(y{R; &) in the polar basis has also 



been determined in ref. [ |130 by a Fourier expansion which is strictly periodic in R. In the first 
period \R\ < L, the series can actually be resummed to give the simple result 



He 



^1 

L 



l = Vc{R)l; 



(3.9) 



for \R\ > Lit must be continued periodically. Quite surprisingly, this turns out to be independent 
of the gauge field t3- and can thus be pulled out of the gauge-fixed expectation value eq. (pj|). 
Inserting again the complete set of eigenstates = of H^^, we arrive at 



(W) =tr(/?|e"^('^fl''+^c-£o)|;j^ 



'tr\{0\U{R 



(3.10) 



The prefactor Zw = tr |(0| ?/(/?) |0)p can be computed explicitly and turns out to T-independent 
(and finite). Again, is the OPE renormalization constant for the composite operator W and we 
obtain the renormalized Wilson loop (Wr) = since Vc{R) is periodically continued for 

\R\ > L, we have to go to the limit L » /? to avoid finite size effects, and the effective string tension 
becomes 



dVc{R) 



dR 



(3.11) 



R=0 



This agrees with eq. ( p3| ) and is also the known expression from the literature. It remains to explain 
why the physical potential Vc{R) is periodic in R, but not invariant under /? — )• L — /? as it should 
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be on a compactified [0,L] ~ 5^ A careful analysis [10] along the lines in sec. ^reveals that the 
physical potential between static quarks is not given by any single Wilson loop once the space 
direction is compactified. Instead, all the equivalent loops with /?—)•/? + mL {m S Z) compete and 
the one giving the minimal potential dominates in the limit T — oo. Taking this prescription into 
account, the level crossing between \R) and \R — L) a.tR = L/2 yields a potential V{R) that is both 
invariant and symmetric about R = L/2. 

3.3 Yang-Mils Theory in D = 3 + 1 



For the physically most intersting case of YM theory in D = 3 + 1 , we start again from eq. ( p.3| ) 
and introduce the usual Coulomb gauge fixing in the Hamiltonian, combined with the resolution of 
GauB' law in the qq sector to which the Wilson state |R) belongs. 



{W) = tr (R_L|e""^^^™+'^'^'''''>'"+''l"'^"^|Rj 



(3.12) 



In view of the approximations below, it is expedient to perform the unitary transformation from sec- 
tion |2| and reshuffle the parallel transporter ?7[A-'-](R) from the Wilson state into the Hamiltonian. 
As a result, we have 



{W) = tr{0\e' 



VHc [Pdyn +Pind +P] ) 



|0) =tr(0|e"^(^fl^"^°)|0), (3.13) 



where the momentum 11 in the gauge-fixed vacuum Hamiltonian is shifted by the induced elec- 



tric field eq. ( |2.5| ), and the Coulomb term eq. ( |2.9| ) will now contain an induced charge pi^d 



[—A-'-, £-'-], in addition to the dynamic charge Pdyn = [—A-'-, 11^] of the gluon and the external 
charge p of the static quarks. To proceed, we decompose the complete Hamiltonian //gx from 
eq. ( |3.13| ) into three pieces, //fix = H^^ + Vc + A//, where H^^ is the usual Coulomb Hamiltonian in 



the absence of external charges, Vc is the non-Abelian Coulomb interaction between the external 
charges, and A// containing the remaining (indirect) charge interactions. 



In order to evaluate eq. ( 3.13 ) with this decomposition, we have resort to a sequence of ap- 
proximations, which are not all under good control: 



1. Jensen's inequality: (W) > trexp 



-r(o| Vc+A//|o) 



-r(0|^'^''-£'o|0) = trexp 

2. Abelization: [U^ (R) , T"] « 

3. Factorization: {0\Q" F"'' Q"'\0) ^ {0\Q" Q"'\0) (0|f "'^|0) = {0\Q"Q'"\0)F 

Approximation (1.) allows to drop H^^, while (2.) removes terms involving the conjugate momen- 
tum n and and (3.) simplifies the remaining pieces in AH and, in particular, the Coulomb 
interaction, Vc- After a lengthy calculation [|^], the Wilson loop is thus reduced to a form 
which requires, as only input from the ground state, the gluon propagator (0|A.'(x)Ay(y)|0) = 
d"'' tij D (x - y) and the vev. of the Coulomb kernel, (0 1 f (x - y ) 1 0) = 5"^ f (x - y ) : 

{W) ssexp. 



T^-Nc I d\x,y)ei{x) Nc'5{x-y)+F{x-y)D{x-y) £,|(y)J>x 



X exp ■ 



TC2g^ 



F{R)-F{0) 



(3.14) 
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The first term in the exponent is again the self-energy of the Wilson lines which gives rise to 
a renormalization Zw of the loop operator. The large distance behaviour of the second term 
is dominated by the small momentum form of the correlators, for which we use the relations 
D{k) ~ [k^ + Myk'-]—2 and F{k) ~ ^ suggested by both variational calculations [Q] and recent 
lattice simulations It can then be shown that the second piece in the exponent of eq. ( |3.14 ) 
is subdominant (~ /?~^) at large distances. The dominant contribution to the Wilson potential is 
given, in our approximation, by the explicit Coulomb interaction, and the Wilson string tension 
equals the Coulomb string tension, Ow ~ Oc- 

4. Conclusions 

In this talk, I have presented a recent calculation of the Wilson potential in the Hamiltonian 
approach to Yang-Mills theory. A careful analysis of GauB' law and the physical degrees of free- 
dom combined with unitary transformations to take care of time evolution reveal that the spurious 
divergences in the Wilson potential can be absorbed in OPE renormalizations of the composite loop 
operator. Using this formalism, the conect result could be derived in cases where the exact ground 
state is known. For the physically most interesting case of Yang -Mills theory in D = (3 + 1), only 
an approximate solution could be obtained in which the Wilson and Coulomb string tension equal 
(while Gw/oc — 0.3 . . . 0.5 in lattice simulations). The main reason for this discrepancy is most 
likely the approximate factorization used in the calculation, which neglects the screening of exter- 
nal charges by gluons. One-loop calculations of this contribution indeed indicate a reduction of 
Ow, and further investigations of this issue are in preparation. 
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